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One-Fe versus Two-Fe Brillouin Zone of Fe-Based Superconductors: 
Creation of the Electron Pockets via Translational Symmetry Breaking 
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We investigate the physical effects of translational symmetry breaking in Fe-based high- 
temperature superconductors due to alternating anion positions. In the representative parent com- 
pounds, including the newly discovered Fe- vacancy-ordered K0.8Fe1.6Se2, an unusual change of or- 
bital character is found across the one-Fe Brillouin zone upon unfolding the first-principles band 
structure and Fermi surfaces, suggesting that covering a larger one-Fe Brillouin zone is necessary in 
experiments. Most significantly, the electron pockets (critical to the magnetism and superconduc- 
tivity) are found only created with the broken symmetry, advocating strongly its full inclusion in 
future studies, particularly on the debated nodal structures of the superconducting order parameter. 

PACS numbers: 74.20.Pq, 74.25.Jb, 74.70.Xa, 75.25.Dk 



One confusing/puzzling aspect of the new iron-based 
high-temperature superconductors is the dilemma of one- 
Fe vs. two-Fe description, concerning the translational 
symmetry of the system. The generic crystal structure 
of these materials consists of two inequivalent Fe atoms, 
distinguished by the alternating tetrahedral coordination 
of the pnictogen or chalcogen anions (c.f.: Fig. la). Since 
this coordination is known to impose dramatic impacts 
on the hopping integrals of Fe d orbitals [1 , the associ- 
ated broken translational symmetry (from 1-Fe perspec- 
tive) is expected to be physically significant and should 
be fully incorporated via the use of the unit cell that con- 
tains explicitly two Fe atoms. Yet, the observed neutron 
scattering intensity [2-5 shows little (if any) indication 
of such broken symmetry; it appears to follow simply 
the 1-Fe Brillouin zone (BZ) of a simple square lattice 
of Fe atoms (Fig. lb). Furthermore, out of convenience, 
most theoretical studies of superconductivity to date do 
not account fully for this broken translational symme- 
try, disregarding the rigorous symmetry constraint. It is 
thus important and timely to clarify quantitatively vari- 
ous aspects of the effects from this broken symmetry (its 
relevance/irrelevance), and to settle, once for all, the con- 
fusing status of the field on the 1-Fe vs. 2-Fe perspective. 

In this letter, three striking effects of the translational 
symmetry breaking potential (TSBP) are revealed by un- 
folding the ab initio electronic band structures (EBSs) 
and Fermi surfaces (FSs) of representative parent com- 
pounds back to the 1-Fe BZ: i) The folded Fe bands 
("shadow" bands) possess overall weak spectral weight, 
explaining the 1-Fe perspective advocated by the neutron 
measurements, and indicating the necessity of the larger 
1-Fe BZ in angle resolved photoemission spectroscopy 
(ARPES) as well, ii) The folding of the bands induces 
an unusual parity switching in their orbital characters, 
suggesting a change of photon polarization in ARPES. 
iii) Most strikingly, the widely discussed electron Fermi 



pockets around (tt, 0) and (0, tt) for supporting super- 
conductivity would not have existed without the TSBP. 
This advocates strongly the full inclusion of TSBP (the 
2-Fe perspective) in theoretical understanding of super- 
conductivity in these materials, and suggests a critical re- 
examination of the debated nodal structure of the super- 
conducting order parameter [6-14 on the electron pock- 
ets. 

Our theoretical analysis is based on unfolded first- 
principles EBSs and FSs of the normal state in 1-Fe 
perspective, which reveals explicitly various aspects of 
the TSBP effects. Standard density functional theory 
(DFT) calculations [15, 16 were conducted with local 
density approximation in the minimum unit cell (8 Fe for 
K0.8Fe1.6Se2 and 2 Fe for the rest). Based on the DFT 
results, symmetry-respecting Wannier functions [TT] with 
Fe d and anion p characters were constructed to capture 
the low energy Hilbert space within [-6, 3] eV, based on 
which the low energy effective tight-binding Hamiltoni- 
ans, H^ were calculated. Finally, unfolded EBSs and FSs 
were obtained via the recently developed first-principles 
unfolding method [18]. 
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FIG. 1: (Color online) Illustration of (a) one-Fe and two- 
Fe unit cells with +(— ) anions located above (below) the Fe 
plane, and (b) the corresponding first Brillouin zones. 
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FIG. 2: (Color online) Unfolded BBSs (top panel) and FSs at /c^ = (middle panel) and kz = tv (bottom panel) in 1-Fe BZ of 
(a) LaFeAsO, (b) BaFe2As2, (c) LiFeAs, (d) FeTe, and (e) K0.8Fe1.6Se2 with spectral intensity colored by green (Fe d^^2_^2, 
d^2_y2^ and dxy)^ blue (Fe dxz)^ red (Fe dyz) and grey (anion p). The dashed lines marks the strict BZ boundaries. The arrows 
illustrate the folding of bands, and dashed ellipse in (c) points to an example of gap openings. 



The basic idea of our unfolding method [TS^ is to simply 
represent the energy-, cj-, dependent one-particle spectral 
function of the real systems (2-Fe zone) using the ba- 
sis from a more symmetric reference system (1-Fe zone): 

^/cn,fen(^) = Y.Kj\{^^\^J)?^KJ,Kj{(^), whcrC K/k 

denotes the crystal momentum of the original/reference 
system, J the band index, and n the Wannier orbital 
index. This change of basis is made simple with the 
use of first-principles Wannier functions [18] . As demon- 
strated below, the unfolded EBSs and FSs provide ex- 
plicit and detailed information on each band's coupling 
to the TSBP in an orbital specific manner. Addition- 
ally, it can be shown [18 that the unfolded spectral func- 
tion corresponds directly to the intensity of ARPES, as 
it includes the main matrix element effects of the mea- 
surement (expect the remaining atomic dipole matrix 
element to be determined based on the chosen photon 
polarization). This use of "regular" momentum distin- 
guishes our method from the glide symmetry-based un- 
folding employed by, for example, Andersen and Boeri 
[19], in which the twisted geometry does not have direct 
correspondence to the ARPES. Similarly, the use of reg- 
ular momentum is essential in the widely applied spin 



fluctuation studies [BHS] of superconductivity via mag- 
netic susceptibility, x(g, cj), since the momentum transfer 
q concerns the difference of two k points. 

Our resulting unfolded EBSs and FSs of the repre- 
sentative families in the nonmagnetic state are shown 
in Fig. 2, colored to emphasize the essential Fe d or- 
bitals. (Enlarged figures focusing around the Fermi en- 
ergy are given in the supplementary materials [16].) A 
few generic features of unfolding can be clearly observed, 
for example, in Fig. 2c. The most obvious one is the 
appearance of the shadow bands, generated from band 
"folding" via the TSBP. Since here the TSBP is of mo- 



mentum q 



TSBP 



= (7r,7r,0) in the 1-Fe BZ unit (except 



for BaFe2As2 and K0.8Fe1.6Se2, whose double layer struc- 
ture gives q^^^^ = (tt, tt, tt) instead) each band is folded 
from k to k -\- q^^^^ ^ as illustrated by the arrows in 
the lower panel. Note that the conservation of spec- 
tral weight dictates a weaker spectral weight for those 
"main" bands that develop stronger shadow bands. Also 
associated with the shadow band formation are the addi- 
tional gap openings occurring at the 2-Fe BZ boundaries, 
indicated by an ellipse in Fig. 2c as an example. Obvi- 
ously, the intensity of the shadow bands and gap opening 



size reflect (although not necessarily represent fully) the 
bands' coupling to the broken symmetry. 

Fig. 2 shows clearly that the anion bands within [-6, - 
2] eV develop very strong shadows bands, of comparable 
intensity to the main bands. This reflects their strong 
coupling to the TSBP, as it is the alternating position- 
ing of the anion that breaks the translational symmetry. 
Given that one can hardly distinguish the main bands 
from the shadow bands, it is obviously more convenient 
to consider these anion bands in the 2-Fe BZ. 

In great contrast, the Fe bands near the Fermi level 
have rather weak shadow bands. In fact, if it weren't for 
the gap openings (some of which are quite large), the Fe 
bands would have looked just like those from a simple 5- 
band system. The overall weak intensity of shadow bands 
explains why neutron spectra appear to respect the 1- 
Fe BZ: Even though the real symmetry of the system 
dictates the 2-Fe BZ, the folding of the spectrum is just 
not strong enough for a clear experimental identiflcation. 
In fact, the lack of folded bands was also reported in a 
recent ARPES experiment. [2^ Consequently, a larger 
1-Fe BZ is necessary in future ARPES measurements, 
since only about half of the EBSs are clearly observable in 
the 2-Fe BZ, where most ARPES to date were conducted. 

Fig. 2 also reveals an interesting orbital-parity switch- 
ing of the band folding, obvious from the change of color 
of the Fe-bands. Consider the FSs, for example. The 
blue (dxz) and red (dyz) bands of odd parity w.r.t. the 
Fe plane always have green {d^z'^-r'^^ dj;2_y2, and dxy) 
shadow bands of even parity, and vice versa. This can be 
understood from the structure of TSBP in these systems. 
Tab. 1 gives the nearest neighbor hopping integrals for 
the low-energy Fe d bands in BaFe2As2 after integrating 
out As p orbitals. It shows that the alternating position- 
ing of the anion (c.f. Fig. la) leads to an alternating sign 
of all teven,oddi ^ud thus breaks the translational sym- 
metry. Consequently, these terms form the main body 
of the TSBP, and dictate a switching of parity in the 
band character upon band folding. This novel behavior 
is quite distinct from the common cases of ARPES, in 
which the replica of bands beyond the flrst BZ retain the 
orbital character. Here, the weak replica always possess 
a different character across the BZ boundaries and thus 
require a different photon polarization to clearly observe, 
similar to the recent reports on Bi2Sr2CaCu208+5 [21 . 

The most signiflcant feature revealed in Fig. 2 is the 
incompleteness of the unfolded electron pockets around 
the X (and R) points. Taking BaFe2As2 in Fig. 2b as 
an example, near the X=(7r,0,0) point, the intensity of 
the strong red pocket vanishes dramatically near the F- 
X path. Consequently, only the green shadow pocket 
of dxy character, folded from the blue pocket around 
R=(0, TT, tt) (c.f. bottom panel of Fig. 2b), is visible here, 
in agreement with recent ARPES measurement [22]. In 
fact, none of the unfolded pockets near X and R points 
are complete in Fig. 2, contrary to the common assump- 



TABLE I: Nearest-neighbor hopping integrals (in eV) along 
the x-direction (r/ + (100), n|if |r/, n/) among Fe d Wannier 
orbitals for nonmagnetic BaFe2As2. The option in sign cor- 
responds to two inequivalent Fe sites. Bold font highlights 
orbitals with odd parity. 
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tion that each X/R point has one strong complete pocket 
and obtains a weaker shadow pocket via band folding (c.f. 
Fig. 3a). All the electron pockets here are instead formed 
by combining main bands near X and R points via the 
TSBP. In other words, without breaking the 1-Fe trans- 
lational symmetry, the essential electron pockets would 
have never existed in these systems. 

To better illustrate this important flnding, let's con- 
struct a translational symmetric Hamiltonian (in 1-Fe 
unit) via the virtual crystal approximation (VCA) to the 
above effective Fe (i-band Hamiltonian: 



H. 



VGA 



i:0,n/ 



/ ^ -tJ-rf-\-r,n;rf,nf/ / ^ I7 



where r denotes the 1-Fe unit cell index. This effectiv- 
ley zeros out all TSBP (the above teven,odd terms), and 
keeps all translational symmetric terms intact. The re- 
sulting EBS and FSs at /c^ = are given in Fig. 3c. Even 
though the overall EBS still follows the original structure 
in Fig. 3b (minus the shadow bands and gap opening ob- 
viously), the topology of the FSs is drastically modifled. 
In particular, there are no electron pockets around the 
X point anymore. The same is found in the VCA of our 
original Hamiltonian containing Fe d and As p orbitals 
(Fig. 3d). Evidently, the TSBP is instrumental in cre- 
ating the electron pockets, and thus should be fully in- 
cluded in future theoretical modeling of magnetism and 
superconductivity. 

Our flndings have direct and signiflcant implications on 
the heatedly debated issue of nodal structures of the su- 
perconducting order parameter on the electron pockets. 
Current spin fluctuation theories [H [HI E] suggest acci- 
dental nodes in the 5^ order parameter on the electron 

pockets (see Fig. 3a), due to strong inter-electron-pocket 
scattering. While the existence of nodes appears to be 
supported by the interpretations of the penetration depth 
and several other measurements [ToHIl], it contradicts 
with the nodeless and almost isotropic gaps observed in 
ARPES [131II5. We flnd that precisely near the region of 
the nodes, ARPES would have negligible intensity, and 
thus can easily miss the nodal structure. On the other 
hand, the above mentioned theories did not incorporated 
appropriately the essential TSBP, and consequently are 
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FIG. 3: (Color online) (a) Folding of complete electron pock- 
ets (from Ref. jGj) and illustration of currently proposed nodal 
structure of the superconducting order parameter, (b) First- 
principles results showing incomplete electron pockets in- 
stead. (c)(d) Demonstration of loss of electron pockets by 
dropping TSBP in both 5- and 8-band VGA descriptions. 



based on FSs of qualitatively different spectral intensity 
and orbital structures. Specifically, one would expect 
that the above incompleteness of the electron pockets, 
the mismatch in the orbital characters, and the addition 
of non-diagonal coupling between the pockets, can all af- 
fect quite strongly the inter-pocket scattering and alter 
the position or even the existence of the calculated ac- 
cidental nodes on the electron pockets. Thus, a careful 
re-examination of the theoretical prediction would be of 
great interest and importance. 

Finally, let's consider the second TSBP introduced by 
ordered 20% Fe vacancies in K0.8Fe1.6Se2 [23 . The long 
period of the vacancy ordering, corresponding to a small 

BZ, making it difficult to compare standard DFT results 
[2T with the experiments. Our unfolded EBS and FSs 
thus offer a direct comparison with experimental spectra. 



and provide a few theoretical insights. Contrary to the 
above discussions, the Fe- vacancy induced TSBP causes 
a strong coherent scattering of the Fe bands, producing 
a larger number of shadow bands and strong gap open- 
ings all over the 1-Fe BZ. Consequently, the overall band 
width of the Fe bands reduces by about 1/4, much more 
than the anion bands. Not surprisingly, the resulting FSs 
are seriously reconstructed, showing little resemblance to 
the generic FSs of other cases, particularly lacking the 
nesting of the Fermi pockets. This disagrees drastically 
with recent ARPES experiment, which reported well- 
defined Fermi pockets [25 . Thus, the Fe- vacancy must 
order only weakly (or locally) in the measured samples, as 
observed recently by electron microscopy [26 . Similarly, 
now it seems more obvious that the recent experimental 
finding of enhanced superconductivity by promoting dis- 
order of Fe- vacancies ^27^ is mostly through the recovery 
of approximately nested Fermi pockets. 

In conclusion, our first-principles unfolded EBSs and 
FSs reveals three key features of the translational sym- 
metry breaking due to alternating anion positioning, in 
all families of the Fe-based superconductors. First, the 
folded shadow bands have rather weak spectral weight. 
This explains the apparent respect to the 1-Fe BZ in neu- 
tron measurements, and highlights the necessity of cov- 
ering the larger 1-Fe BZ in future ARPES experiments. 
Second, TSBP induced band folding changes the orbital 
character to those with opposite parity w.r.t. the Fe 
plane. This unusual phenomenon suggests a change of 
photon polarization in ARPES experiment. Finally and 
most significantly, the electron pockets, critical to most 
theories of magnetism and superconductivity of these ma- 
terials, only form via the TSBP. Thus, full inclusion of 
the broken translational symmetry (e.g. using 2-Fe unit 
cell) is essential in future theories, particularly on the 
debated issue of nodal structure of the superconducting 
order parameter on the electron pockets. 

Work funded by the U S Department of Energy, Office 
of Basic Energy Sciences DE-AC02-98CH10886 and by 
DOE-CMCSN. 
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DETAILS OF AB INITIO CALCULATION AND BAND STRUCTURES 

All ab initio calculations were conducted by density functional theory with local density approximation in linearized 
augmented plane wave basis, which is implemented in WIEN2k package [1]. For each parent compound, the lattice 
constants and atomic positions were obtained from the experimental data of the high temperature nonmagnetic state of 
LaFeAsO [2 , BaFe2As2 [3 , LiFeAs [4 , FeTe [5^, and K0.8Fe1.6Se2 [6 . The essential parameters used in the calculation 



are summarized in Tab. 2. We followed the default settings of version 10.1 with RmtKn 



7 and L 



10 to 



reach convergence of the ground state density. Then the symmetry respecting Wannier functions [11 18] ofFed and As 
p orbitals were constructed within the low energy Hilbert space from -6 to 3 eV. Thus, the orbital-resolved spectral 
function in the 1-Fe Brillouin zone [9] can be calculated from the effective tight-binding Hamiltonian on Wannier 
orbital basis by the unfolding method \9 . Especially for the two-dimensional Fermi surfaces in kz= and kz = tt 
planes in the letter, we plotted the spectral function on the 400 x 400 k point grid in order to capture all fine features. 
The unfolded band structures focusing on ±2 eV energy around Fermi energy were plotted in Fig. 4 to reveal the 
details of Fe d bands around Fermi energy. 



DETAILS OF VIRTUAL CRYSTAL APPROXIMATION 



In the presence symmetry breaking potential, 
virtual crystal approximation (VGA): 



a translationally symmetric Hamiltonian can be constructed via the 



TtVCA _ \~^ TT 

-"-r,n;0,n/ — / ^^r 



f-\-r,n]rf,nf / / ^ I7 



(1) 



where Hru^r.n^n.ni and HY^-Qn/ ^^^ Hamiltonians with broken translational symmetry in the basis of Wannier orbitals 
n located at normal cell index r. The summation of r/ is over inequivalent normal cell lattice vectors in a single 
supercell (namely the positions of two or more inequivalent Fe sites). Thus, the symmetry breaking terms of the 
original Hamiltonian would be averaged out, but the symmetry respecting terms would be kept. 

In the Fe-based superconductors, one can choose whether to integrate the As p orbitals into Fe d Wannier orbitals 
or not, during the construction of the Wannier orbitals. If As p orbitals are integrated out, we will effectively obtain 
the five-band VGA Hamiltonian from Eq. 1. If not, the eight-band VGA Hamiltonian (5 Fe d and 3 As p) will be 
obtained. Fig. 3 in the letter shows clearly that VGA has a larger impacts in the 8-band case. This is easily understood 
from the following consideration. In the 8-band case, the symmetry breaking terms are mostly those hopping terms 
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ABLE U: Lattice constants and k meshes used 


in the density functi 


onal theory calculation. 




Lattice constants 


LaFeAsO 


BaFe2As2 


LiFeAs 


FeTe 


K0.8Fe1.6Se2 


a=b (Bohr) 
c (Bohr) 


7.615732 
16.51017 


7.488039 
24.59819 


7.127100 
12.01380 


7.221591 
11.84764 


16.47577 
26.61737 


k mesh 


16 X 16 X 7 


13 X 13 X 13 


14 X 14 X 8 


17x 17x 10 


11 X 11 X 11 
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FIG. 4: Unfolded electronic band structure in 1-Fe Brillouin zone of (a) LaFeAsO, (b) BaFe2As2, (c) LiFeAs, (d) FeTe, and 

(e) K0.8Fe1.6Se2 with spectral intensity colored by green (Fe d^z^_r2, dx2_y2, and dxy), blue (Fe dxz), red (Fe dyz), and grey 
(anion p) . 

involving Fe d and As p orbitals, since As p orbitals are located in different locations. Upon integrating out these 
terms, they effectively renormalize the remaining hopping terms between Fe d orbitals, both the symmetry breaking 
and symmetry respecting ones. By throwing out the symmetry breaking terms at the eight-band level via VGA, one 
also removes the renormalization to the symmetry respecting hopping between Fe d orbitals as well. In either case, 
however, the qualitative effect of losing the electron pockets persists, since it is mostly dictated by the symmetry 
property of the system and is thus more robust against detail changes of the Hamiltonian. 
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